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1 A simple two-layer model, the moist-convective rotating shallow water, which al- 
2 lows for low-cost high-resolution numerical simulations of the dynamics of the moist 
3 atmosphere in the presence of topography, is used to identify and understand dy- 
4 namical processes governing the evolution of easterly waves propagating on the 
5 background of a low-latitude easterly jet crossing a land-sea boundary, a setup 
6 crudely representing the African Easterly Jet over the West-African plateau and 
7 the Atlantic ocean. We perform a thorough linear stability analysis and identify the 
8 unstable modes of the jet, which we use then for initialisation of fully nonlinear 
9 numerical simulations. In this way we determine nonlinear evolution of unstable 
10 perturbations of the jet, both in the “dry” and moist-convective environments, and 
i highlight essential differences between the two cases. We identify a mechanism of 
12 formation of intense lower-layer cyclonic vortices at the northern flank of the jet, 
13 and determine the influence of the land-sea contrast upon this process. 

14 Keywords: Baroclinic instability, Moist convection, Orographic effects, African Easterly Jet, Tropical 

15 Cyclogenesis 
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1. Introduction 


Understanding tropical cyclogenesis in the North Atlantic is of obvious importance. It is 
well-known that the cyclogenesis is most often triggered by depressions due to the African 
Easterly Waves (AEW) propagating over Atlantic, (e.g. Bracken and Bosart 2000, Brammer 
and Thorncroft 2015, Dieng et al. 2017), and references therein. In turn, it is widely accepted 
that AEW are related to instabilities of the African Easterly Jet (AEJ), a low-latitude atmo- 
spheric easterly jet evolving over the West-African plateau and the Eastern Atlantic Ocean 
(Burpee 1972, Hsie and Cook 2005, 2008, Wu et al. 2012). Understanding the instabilities of 
AEJ is, thus, of primary importance, and their studies, which started almost half a century 
ago (Simmons 1977) still continue (Thorncroft and Hoskins 1994a, Hall et al. 2006, Diaz and 
Aiyyer 2015). It is known that humidity, and the contrast between the ocean and the continent 
play an important role in the development and intensification of AEW (e.g. Brammer and 
Thorncroft 2015). 

In the present paper, we analyze the development and nonlinear evolution of the insta- 
bilities of low-latitude easterly jets in the presence of moist convection and the land - sea 
contrast with the help of a simple two-layer atmospheric model which, nevertheless, captures 
essential dynamical mechanisms. An advantage of the model is that inclusion of topography 
is simple and straightforward. The model allows for efficient low-cost high-resolution linear 
and nonlinear analyses of the system, and we profit from this to establish fine details of lin- 
ear stability diagrams, and of nonlinear saturation of instabilities. An early precursor of this 
study is the work by Kuo (1978), where barotropic and baroclinic instabilities of easterly 
jets were analyzed, with similar motivations, in the framework of a two-layer model of the 
same kind, although only on the linear level and with inclusion of neither topography nor 
dynamical effects of moisture and condensation. The model we are working with is the so- 
called moist-convective rotating shallow water (mcRSW), which, in its two-layer version, was 
derived in Lambaerts et al. (2011b), by using vertical averaging between isobaric surfaces and 
Lagrangian conservation of the equivalent potential temperature in the hydrostatic primitive 
equations with pseudo-height as vertical coordinate. This conceptual model proved to be use- 
ful in studies of instabilities of large- scale Earth and planetary jets and vortices, and the effect 
of diabatic effects upon them (Lahaye and Zeitlin 2016, Rostami and Zeitlin 2017, Rostami 
et al. 2017, 2018). The main goal of the present paper is to show that, in spite of a highly 
idealized character of the model, it allows to capture the fundamental dynamical mechanisms 


governing the evolution of AEW. Our strategy below is similar to that of Lambaerts et al. 
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4 M. Rostami and V. Zeitlin 


(2012), where the model, for benchmarking purposes, was applied to the moist baroclinic in- 


stability of the upper-tropospheric midlatitude westerly jet on the f- plane. Yet the important 





differences, apart from the jet configuration itself, are the use of the beta-plane approximation 
and addition of topography with a land-sea contrast. We choose a basic state which mimics 
the AEJ, within the precision of the model, and perform, first, the linear stability analysis 


in the “dry” case, with condensation switched off, and without topography. In this way we 





identify the unstable normal modes of the jet, which we use then, either in a pure harmonic 
form, or in a form of localized wave-packets, to initialize fully nonlinear high-resolution nu- 
merical simulations, which are performed both in “dry” and moist-convective environments, 
and also include an idealized step-like topography mimicking the West-African plateau and 


the land - sea contrast. We are able in this way to highlight the role the moist convection 





and topography are playing in intensification of the low-level depressions created by unstable 
disturbances, and their eventual transformation into strong localized vortices. 

The paper is organized as follows. In section 2 we present the model, introduce an appro- 
priate scaling, and define the background jet configuration. In section 3 we give the results 
of the linear stability analysis of the jet configuration in the “dry” version of the model, with 
the condensation switched off. In section 4 we present results of numerical simulations of non- 
linear evolution of the unstable disturbances in the “dry” and moist-convective environments, 
both without and with topography. Section 5 contains summary and discussion. We give in 
Appendix A the results of the stability analysis in the f-plane approximation, which was used 


for benchmarking purposes. 
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INSTABILITIES OF EASTERLY JETS 5 


2. The model, the scaling, and the jet configuration 


2.1. The moist-convective rotating shallow water (mcRSW) model 


We will work with the later version of the mcRSW model (Rostami and Zeitlin 2018) defined 


by the following set of equations: 


y-7V v1 v2 


Ov, + (01° V)ur t+ fly)Z x vy = -—gV (hi + he) + 7 ; (2.1a) 
Opv2 + (vo: V)va + fly)Z x ve = —gV (hi + sho) + eee SS (2.1b) 
Ophy + V - ((hy — b)v1) = —yCO +-7°V, (2.1c) 
Oho + V - (hove) = +7C — y°V, (2.1d) 
0.Q, + V.(Qiv1) = -C+E, (2.1e) 
0:,Q2 + V.(Qov2) = V. (2.1f) 
OW, + V.(W101) = +(1 -—«)C- P, (2.1¢g) 
O,We + V.(Wove) = +k 0 —V, (2.1h) 


where v; = (u;, v;) is the horizontal velocity field in the layers i = 1,2 counted from the 
bottom, f(y) = fo + By is the Coriolis parameter, Z is a unit vector in the vertical direction, 
h; are the thicknesses of the layers in the absence of topography, b represents the topography, 
g is the gravity acceleration, and s = 62/0; > 1 is the stratification parameter, where 6; are 
constant potential temperatures of the layers. Without the underlined terms (2.1) becomes the 
standard non-dissipative two-layer atmospheric shallow-water model (cf. Zeitlin 2018) with an 
addition of topography and two passive scalars layerwise: the specific humidity (water vapor 
content) integrated over the air column Q;, and the bulk amount of liquid water W; in the 
column. Both are conserved quantities if diabatic effects are switched off. We will call the 
system “dry” in the latter case. The underlined terms introduce a condensation sink C’ and 
a surface evaporation source F in the water vapor equation in the lower layer (2.le), and a 
convective flux between the layers in the mass conservation equations (2.1c), (2.1d), due to 
the latent heat release. ‘his flux is proportional to the condensation, with the proportionality 
coefficient y which is defined by the underlying vertical stratification and thermodynamical 
parameters (cf. Lambaerts et al. 2011b). Similar parameterization is used for vaporization V, 
which is a source of water vapor, and introduces cooling with a related downward convective 


flux which is proportional to vaporization, with a proportionality coefficient y*. We suppose 
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6 M. Rostami and V. Zeitlin 


that the lower layer is close to saturation, and neglect vaporization in it and that the upper 
layer is far from the saturation and neglect condensation in it. Correspondingly, condensation 
(vaporization) is a source (sink) of liquid water. Entrainment of liquid water from lower to 
upper layer due to the upward convective flux is also introduced, and is controlled by the 
parameter «. Precipitation P is a sink of liquid water. The upward and downward convective 
fluxes lead to the Rayleigh drag terms appearing in the right-hand sides of (2.1a), (2.1b). 

Relaxational parameterizations are used for the condensation C' and precipitation P, with 
corresponding relaxation times T., Tp: 

W — Wer 


_&-“yW9g_g), pa —"* ww _w.), (2.2) 


Te Tp 


C 


where Q* is the saturation threshold for water vapor, W-, is the precipitation threshold, and 
H denotes the Heaviside (step-) function. Bulk formulas are used for surface evaporation and 


vaporization: 


v1 | (Q° — Q) 
Qs 


where |Umaz| is the maximum value of the lower-layer velocity, and ae, are adjustable coeffi- 


B=-a 





(Q° —-Q)H(Q*-Q), Vay WH(Q* — Q) (2.3) 


e 
Wniax | 


cients. Although a threshold for evaporation can be taken to be different, we make the simplest 


choice and set it to be the same as in the condensation (2.2). The saturation threshold Q* 





depends on pressure according to the Clausius-Clapeyron law. This dependence, within the 
approximations of the model, can be taken to be linear (cf. Bouchut et al. 2009) or exponen- 


tial, with a small exponent: Q*(h) = Qje  & Q5(1 — en), where 7 is the thickness anomaly 





in the lower layer, see below. Evaporation (2.3) is the only process in the boundary layer we 
include in the model. Momentum dissipation due to the bottom friction can be easily included, 
and more involved parameterizations of the boundary layer can be used too (cf. Schecter and 
Dunkerton 2009), which we will not do, keeping the model as simple as possible and limiting 
the number of adjustable parameters. The parametrization of surface evaporation in (2.3) is 
well suited for the evolution of the atmosphere over the ocean (e.g. Katsaros 2001) but should 
be changed over land. Again, for simplicity, we simply switch it off whenever b 4 0. We should 
emphasize that an advantage of this version of the model, as compared to the original one 
in Lambaerts et al. (2011b), is that it tracks condensed water, and thus is “cloud-resolving” , 
with the distribution of W giving a rough idea of the cloud cover layer-wise. 

Radiative relaxation term of the form —(h; — H,)/7;, where 7, is a characteristic relaxation 
time, can be added in the r.h.s. of (2.1c) in the condensation regions, to model the physical 


process of enhanced outgoing radiation by clouds. We test the influence of this term in the 
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INSTABILITIES OF EASTERLY JETS 7 


numerical simulations below. 
A detailed discussion of the conservation laws and energy budget in the model can be found 


in Rostami and Zeitlin (2018) and we will not repeat it here. 


2.2. Scaling of the “dry” version of the model, and the background jet configuration 


We introduce unperturbed thicknesses of the layers Ho;, 2 = 1,2 and the unperturbed total 
thickness Hp = Ho; + Hog. The following scaling is used for the independent and dependent 


variables in the “dry” version of the model: 


(x, y) ~la=v gHo/ fo, EO (1/ fo), (u, v) ~ Vg, (h, b) ~ Ho, (2.4) 


where Lg denotes the barotropic deformation radius. For dimensional estimates we will place 
the beta - plane at a low latitude ®p, and take Hp = 10km, fp = 22 sin ®p = 3.2717x10-° ss! 
at &y = 13° N, which gives Lg = 9568km, 3 = 2.2243 x 107" (ms)~!, and /gHo = 313 ms™". 


We recall that in the tangent-plane approximation 


20 L 7 
f = 202 sin@®@ = 222 sin &p + Fb = fo+ Py = fo (: = —ty") = fo(14+ By"), (2.5) 
. 0 


where the asterisk denotes non-dimensional variables, and we introduced the non-dimensional 
beta 8. Its value is 8 = 6.505 at ® = 13°. The non-dimensional “dry” equations on the 


6-plane, after dropping the asterisks become : 


OU; + UjzOr Uj + Vj Oy Ui — (1 ++ By)v; = —O;(h4 + s' tha), (2.6a) 
OnV; + UjO,v; + ViOyv; + (1 + By)u; = —Oy(hi + s' tho), (2.6b) 
Athi + Ory (hg — (2 — 6)b)uy) + Oy ((Ag — (2 — )B)vi) = 0, §= 1,2. (2.66) 


It is easy to see that in the absence of nontrivial topography, i.e. with 6 = const, any zonal 


jet u; = U;(y), hi = Ho; + Ai(y), vi = 0 in geostrophic equilibrium: 
(1+ By)U; = —0,(Hy + s’~* He) (2.7) 


is an exact stationary solution of (2.6). If b = b(2) 4 0, which will be our choice below, only 
the upper-layer jets with quiescent lower layer U; = 0 are exact solutions. ‘This will be the 
case we consider. We choose the classical Bickley jet profile for U2, which we orient from East 


to West: 


U2(y) = —Up sech* [aly — yc)], (2.8) 
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Figure 1. Meridional profiles of zonal velocities U; (left panel) and thickness deviations H;(y) (right panel). Lower 
(upper) layer: solid (dashed). Chebyshev nodes, which are used for spatial discretization, are shown in the right panel. 
Ro = 0.15, Bu = 11.35, Up = 0.045, a = 20, Ho1 = 1/3, Ho2 = 2/3, k = 20, 62/0; = 1.1, 6 = 13°N. 


136 Where a controls the width of the jet. ‘The Burger number of the jet is defined as Bu = L* /R?, 


i37 Where RF is the width of the jet at the level 0.01 Up. ‘The position of the jet axis y, is variable. 





i33 Instead of changing it, we choose y, = 0 and vary instead the latitude @g, and hence the value 
39 of 8 in (2.6), as it is known that this parameter is crucial for (in)stability properties (Kuo 
140 1978). The corresponding profiles of background thickness deviations H; are recovered from 


141 the geostrophic balance 


Hay) =— | eal. Tela): (2.9) 


ao 
v2 An upper-layer jet solution, with the characteristics roughly corresponding to the AEJ (e.g. 
uz Hsie and Cook 2008) is presented in Figure 1. 


4 3. Linear stability analysis of the upper-layer easterly jet in the absence of 


145 topography 


We introduce small perturbations of all variables, denoted by primes, with respect to the 
upper-layer jet solution (2.8), (2.9): 
W= uy, w=v,, k= Ai(y)+7+ Ao, 


ug = Uo(y)+ug, v2= vg, he = Holy) ++ Hoo, 


and linearize the equations (2.6) with b = 0 about the stationary solution. Dropping the 
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primes, we thus get: 


Oyu — (1+ By)vi + O2(m + 2) = 0, 

Oyu + (1+ By)ur + Oy(m + 2) = 0, 

On, + HyOzu1 + v1 0, + Hy0yv1 = 0, 

Oyug + Ug0,ug + V2DyUq — (1 + By)v2 + Ox(m + 8N2) = 0, 
Oyvg + Uzdxv2 + (1+ By)ua + Oy(m + 872) = 0, 


O1n2 + U20zN2 + H20;uU2 + v20yHe + H20yv2 = 0. 


We look for harmonic solutions of (3.10) in the form: (uj;,v;,7) = Re [(%, it, 7) ef 








(3.10a) 
(3.10b) 
(3.10c) 
(3.10d) 
(3.10e) 
(3.10f) 


ay) 


Complex eigenfrequencies w = wr + iw; with positive imaginary part (wy; > 0), correspond to 


instabilities with linear growth rate 0 = wy. ‘The resulting eigenproblem is: 


1 jt V1 inl U2 V2 ro] = w | ti V1 1 U2 V2 me] 


0 —(1+ By) k 0 0 
—(1+ By) 0 —D 0 0 
; kHy Hi(y)+Hi(y)D 0 0 0 
7 0 0 k  —— kU2(y) ~—-U3(y) — (1 + By) 
0 0 =D =(1-By) kU2(y) 
0 0 0 kie(y) Aly) + Ha(y)D 


(3.11) 
k; 
—D 
0 
ks | 
—sD 
kU2(y) 
Cae) 


where D denotes the operator of differentiation with respect to y, which will become the 


Chebyshev differentiation matrix after discretization, the superscript I’ denotes the transpose, 


and prime denotes the derivative of the corresponding functions with respect to their argument 


Y. 


We should emphasize at this point that the analysis of the spectrum of unstable modes 


of easterly jets on the beta-plane, as compared to their westerly counterparts, encounters 





difficulties related to the presence of critical levels y = y,; where the phase velocity of the 


eigenmode c = wp/k coincides with the jet velocity U(y-;), which appear already in the sim- 


plest case of the barotropic jet in the quasi-geostrophic approximation (e.g. Maslowe 1991, 


Swaters 1999). This is why a theoretical analysis of this spectrum is a challenge. Instead of 
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Figure 2. Stability diagrams of the upper-layer easterly Bickley jet on the 6-plane at three different latitudes in terms 
of phase velocity cp (left panel) and growth rate o (right panel) as functions of zonal wavenumber k. Ro = 0.15, Bu = 
11.35, Uo = 0.045 | foLal, a = 20, Ho1 = 1/3, Ho2 = 2/3 (Colour online). 


tempting it, we adopt a pragmatic approach, consisting in a direct numerical solution of the 
eigenproblem (3.11), (3.12) by using pseudo-spectral collocation method with a non-uniform 
distribution of Chebyshev nodes which is denser inside the jet and is shown in the right panel 
of Figure 1. Critical levels should be treated carefully, as they can give rise to discontinu- 
ous numerical pseudo-modes. We identify them by increasing the numerical resolution, and 
discard. ‘Typically, the resolution with N = 400 collocation points in y turns out to be fully 


sufficient. We first benchmarked the method on the f-plane, where the results are insensitive 





to the orientation of the jet on the plane, and had to reproduce those of Lambaerts et al. 
(2012) who studied stability of the same configuration, but in the case of a westerly jet. This 
is indeed the case, both for the stability diagram, and for the phase-portrait and meridional 
cross-section of the most unstable mode, as shown in the Appendix. 

After having benchmarked the method, we perform the same analysis on the beta-plane. 
Stability diagram in a large range of zonal wavenumbers for three latitudinal positions of the 
jet, with corresponding changes of 2, is presented in Figure 2, where we display not only 
the first, most unstable, but also the second unstable mode. As follows from the Figure, the 
higher the latitude the closer, qualitatively, the results are to those on the f-plane, cf. Figure 
Al, although the growth rates on the beta-plane are higher. The curves of the phase velocity 
display a non-monotonicity in the interval k € [5, 20], which becomes more pronounced at 
lower latitudes. A similar non-monotonicity was reported in Kuo (1978), although that paper 
used a two-layer model with a rigid lid and explored a configuration with nonzero jet velocity 
in the lower layer, but also in Thorncroft and Hoskins (1994a), who worked with the primitive 


equation model on the sphere, like Simmons (1977). We should emphasize a very good qual- 
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Figure 3. Real part of eigenfrequency (left panel) and group velocity (right panel) of unstable modes of Figure 2 (Colour 
online). 


itative agreement of the stability diagram, as well as a good quantitative agreement of the 
zonal wavenumber, phase velocity, and growth rate of the most unstable mode in our simple 
model (taking into account our choice of scaling), with the results of Thorncroft and Hoskins 
(1994a). Note that only the most unstable mode was identified in this paper, while we find 


them all. What is important, the corresponding group velocity of the unstable modes, which 





is presented together with the real part of the eigenfrequency in Figure 3, is always westward, 
as the phase velocity itself, like in the case of AEW in the observations. 


The phase portraits of the first and the second unstable modes of Figure 2, with the same 





zonal wavenumber k& = 20, are presented in Figure 4. First, we clearly see in the Figure that 


the phase portrait of the most unstable mode is qualitatively similar to that of its counterpart 





on the f-plane, cf. Figure A2, although it is subject to deformations due to the beta-effect, 
especially in the lower layer. Second, we see that the first unstable mode is varicose, 1. e. 


having anomalies of the same sign at both sides of the jet at a given longitude, and the 





second unstable mode is sinuous in the upper layer, where the jet is located, a general fact 
which is well-known for jet instabilities. Note that baroclinic character of the unstable mode is 


more pronounced at the southern side of the jet. The instability in question can be, as usual, 





interpreted as a resonance between Rossby waves counter-propagating across the meridional 
gradients of potential vorticity (cf. Thorncroft and Hoskins 1994a). An important question 
raised in the literature (Hall et al. 2006) is sensitivity of the results of the stability analysis to 
the parameters of the jet. Within our simple model, and with our simplest parameterization of 


the jet, a study of sensitivity is straightforward. The sensitivity of the results to the latitudinal 





position of the jet is clear from Figure 2: the lower the latitude of the jet center, the larger the 
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Figure 4. Geostrophic stream-functions of the unstable modes on the {-plane in the upper layer (upper row) [w2 = 
71 + s7y2] and the lower layer (lower row) [wi = 71 + 72] of the first (left column) and the second (right column) unstable 
modes. Positive (negative) values are represented by solid (dashed) lines, and the velocity field is represented by arrows. 
Ro = 0.15, Up = 0.045, Bu = 11.35, a = 20, Ho1 = 1/3, Hoo = 2/3, k = 20, ® = 13° N (Colour online). 





wavenumber of the most unstable mode. The sensitivity to the value of the Burger number, 
i.e. to the width of the jet, can be inferred already from Figure Al. The sensitivity to the 
Burger number on the (§- plane is presented in Figure 5. As follows from these figures, at lower 
Bu, i.e. for wider jets, the wavenumber of the most unstable mode diminishes, and at higher 


latitudes the wavenumber of the most unstable mode shifts down, as well. 





The results of the analysis of sensitivity to the aspect ratio and stratification are presented, 
respectively, in Figures 6 and 7, which show that these sensitivities are rather low in the 
interval of wavenumbers corresponding to the highest (and, thus, most interesting) growth 
rates. 

We should finally comment on the role of the bottom friction, which is easy to include in 
the model, as already mentioned. We, however do not do this, as the lower layer in our basic 
state is quiescent, and the amplitude of velocity of the unstable modes in this layer is an 


order of magnitude lower than in the upper layer, cf. Figure 4. Yet in the configuration with 





lower-layer westerlies, which we neglect, the lower layer damping does influence the structure 
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Figure 5. Sensitivity of the stability diagram, in the region of the highest growth rates to the value of Bu. Phase velocity 
Cp (left panel) and growth rate o (right panel). Ro = 0.15, Up = 0.0447 [foLa],k = 20, s = 1.1, Ho1 = 1/3, Ho2 = 2/3. 
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Figure 6. Sensitivity to the ratio of the layers’ thicknesses of the stability diagram for phase velocity cp (left panel) 
and growth rate o (right panel). Ro = 0.15, Bu = 11.35, Up = 0.045 [fo La], a = 20,k = 20, s = 1.1, Ho1 + Ho2 = 1. 


of the normal modes (Hall et al. 2006). 


4. Nonlinear evolution of the unstable disturbances 


4.1. Parameters and setup 


To investigate the nonlinear evolution of the instabilities of the jet, we used a high-resolution 
finite-volume code by Bouchut and Zeitlin (2010). The simulations were set in the domain 


with the zonal dimension L, = nA, n = 5, in units of Lg, where A is the wavelength of an 


September 15, 2021 


224 


225 


226 


227 


228 


229 


230 


251 


232 


233 


234 


255 


236 


237 


238 


239 


240 


241 


Geophysical and Astrophysical Fluid Dynamics 


























LLEJ GAFD'RZ2021' Jun30 


























14 M. Rostami and V. Zeitlin 
-0.03 + 0.16 
-0.031 + 0.14 
-0.0382 F 
0.12 
-0.033 + 
0.1 
-0.034 
om S 
-0.035 0.08 
-0.036 eal 
-0.037 aes 
-m-0,/0, =1.1,1 $* Mode 0.04 -f§——°=e4 -m-0,/0, =1.1,1 St Mode 
-0.038 -O-6,/0, = 1.5,1 °* Mode Lice -O- 6,/0, = 1.5, 1 ** Mode 
BS 
eeeae | = 6,/0, =1.1,2 "* Mode O:027- ~ = 0,/0, =1.1,2 "* Mode 
—o-6,/0, = 1.5, 2 "* Mode Zo |-o 6,/0, =1.5,2 "* Mode 
-0.04 l l l l l 0 l l l ! l 
15 20 20 30 35 15 20 25 30 35 
kK kK 
Figure 7. Sensitivity to stratification of the stability diagram for phase velocity cp (left panel) and growth rate o (right 
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unstable mode used in the simulation, and meridional dimension L, = 3.8L. A “continent” 





in the simulations with topography was represented by a meridional band of width 0.7Lg, 
of uniform height b = 0.05Hpo. (The results presented below do not qualitatively vary with 
changes of this value, within the same order of magnitude). Periodic and sponge boundary 
conditions were implemented, respectively, in zonal (along-jet), and meridional (cross-flow) 
directions. Notice that the periodic boundary conditions imply that the background jet has 
uniform intensity across the domain, which is, obviously, an idealization of the observed AEJ. 
We should also say that, concerning the land-sea boundary, the model, by construction, cannot 
capture the associated thermal contrast. To do this, a variant of the model, which was recently 
proposed in Liu et al. (2020), and allows for horizontal temperature gradients, is necessary. 
This development is postponed to a future work. 

The following values of parameters characterizing the moist processes were used in the 
simulations: y = 0.3, y* = 0.ly, ae = 0.25,€ = 0.3, Q° = 0.9, W., = 0.01. The simulations 
in the moist-convective environment were initialized with a uniform, humidity in the lower 
layer, close to saturation, according to the hypotheses of the model: Q; = Q* — 0.01 over 
the ocean, and far from the saturation 0.2 Q° over the land. ‘The humidity in the upper layer 


was taken to be far from the saturation (Q2 = Q* — 0.2) The liquid water content was set to 








be zero everywhere at the initial moment. ‘The difference between the evaporation properties 
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between the land and the sea was accounted for in the simplest way: the evaporation was 
switched off over topography. The condensation and precipitation relaxation times were taken 
to be T, = 1200; & 2h, Tp = 300:, respectively, where dt = 0.002 is the non-dimensional 
time-step of the code. Increasing the relaxation time of condensation within the same order 
of magnitude does not change the results qualitatively. ‘The non-dimensional vaporization 
coefficient was chosen as a, = (10- @Q*r7,.)~! We also compared simulations without and with 


radiative cooling, with the relaxation time 40 days in the latter case. 


4.2. “Dry” vs moist-convective simulations without topography initialized with the 


most unstable mode 


We start with the analysis of nonlinear evolution of the most unstable mode in the absence 


9 


of topography in the “dry” configuration, by initializing the simulations with the background 
jet with superimposed small (® one tenth of percent, compared to the values in the jet 
itself) perturbation of all fields of an unstable mode. The corresponding evolution of relative 
vorticity is presented in Figure 8, and reveals no characteristic comma-shaped pattern of 


a breaking unstable wave, which is typical for midlatitude westerly jets (Lambaerts et al. 





2011a). Nevertheless, a steepening of the streamlines, which was called “breaking AEW” in 
the analysis of observations (Cornforth et al. 2017, Lafore et al. 2017), does take place. After 
a period of initial growth a nonlinear saturation at finite amplitude takes place. At the same 
time, lateral spreading of vorticity anomalies in the lower layer is observed. What is important 
in the following, this spreading is meridionally asymmetric, extending farther at the northern 
than at the southern flank of the jet. As is easy to check by using the values of the phase 
velocity and frequency from Figures 2, 3, neither of them matches their counterparts for 
free Rossby waves estimated from the quasi-geostrophic expression w = —Bk/(k? + 1? + ie ), 
where / is meridional wavenumber of the free wave, and Lg is deformation radius (barotropic or 
baroclinic). Analysis of the evolution of divergence (not shown) reveals no significant emission 
of inertia-gravity waves either. Hence, the origin of the saturation of the instability is not a loss 
of energy by radiation, and is essentially nonlinear. It is not a breaking of the unstable mode, 
like in westerly mid-latitude jets, either, but a repartition of energy through the spectrum of 
unstable modes. This is confirmed by the Fourier analysis of the dynamical variables, which 
reveals the appearance and growth of modes with lower values of k during the process of 
saturation, not shown. (Notice that subharmonic modes were excluded in Lambaerts et al. 


(2012) by the choice of boundary conditions, unlike the present case). 





As compared to the “dry” simulation of Figure 8, the evolution of the instability in the moist- 
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Figure 8. Evolution of vorticity and velocity in the upper (upper row) and lower (lower row) layer s in “dry” environment 
in a simulation initialized with the full periodic unstable mode of small non-dimensional amplitude (0.005), without 
radiative cooling. The global configuration is the same as in Figure 4. The length unit at the axes here and below is Lg 
(Colour online). 





convective (MC) environment is totally different, as follows from Figure 9. The meridionally 
asymmetric periodic zones of cyclonic vorticity developing from their counterparts in the 
initial perturbation, like in the dry case, undergo a different evolution under the influence of 
the moist convection. At earlier stages they are intensified due to the moist convection. ‘The 
mechanism of this intensification in the present model is very simple, and consists in growth 
of cyclonic vorticity due to the convective mass sink, as explained in Lambaerts et al. (201 1a). 
As follows from Figure 9, the condensation is correlated with the enhanced convergence zone 
inside of the vorticity anomaly, which is shifted north-east. Convergence is a proxy for ascent 
in shallow-water models, so what is observed in the Figure is consistent with the discussion 
of the position of maximum ascent in Thorncroft and Hoskins (1994a,b). At later stages a 
secondary lower-layer jet, which is retrograde with respect to the primary one, is formed 


close to the axis of the original upper-layer jet. Interestingly, retrograde jets in the lower 
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troposphere are typically accompanying the AEJ in the observations (Cook 1999), so here 
we see a dynamical mechanism of their generation. As follows from the Figure, this jet is 
formed by merging cyclonic anomalies of relative vorticity, which are induced by the unstable 
mode, are growing due to the instability, and are enhanced by the condensation. On the other 
hand, localized zones of enhanced cyclonic vorticity are formed at the northern flank of the 
jet. No such cyclogenesis is observed at the southern flank, due to meridional asymmetry 
of the unstable mode. As follows from Figure 10, where we present the late stages of the 
simulation of Figure 9, the cyclones that form at the northern flank of the jet are further 
intensified due to the upward convective fluxes, as explained in Lambaerts et al. (2011a), and 
have a tendency to detach from the jet. ‘Transformation of positive vorticity into detached 
cyclones is not fast and takes few weeks to occur. This suggests an importance of non-modal 
growth, and corresponding initializations of non-linear simulations, see below. Such properties 


of cyclogenesis by AEJ, and asymmetry between the northern and southern flanks of the AEJ, 





have been reported in data analyses (e.g. Pytharoulis and Thorncroft 1999, Chen et al. 2008). 
The cloud cover presented in the bottom-right panel of the Figure is in qualitative agreement 
with observations, in a sense that condensation is within and eastward of the trough, but more 
complicated cloud structures were also detected in the observations, e.g. Cornforth et al. (2017) 
and references therein. 

The left and the right panels of Figure 11 display Hovmoller diagrams of vorticity corre- 


’ and moist-convective experiments, and clearly show 


sponding to the above-described “dry’ 
the westward propagation of disturbances, and their intensification in the moist-convective 
case. 


We also repeated the same simulation with radiative relaxation. Qualitatively, the scenario 





of the evolution does not change, but the intensity of the vortex and jet structures generated 


by the developing instability decreases at later stages, not shown. 


4.3. Evolution of localized wave-packets of unstable modes, and the role of the 


land-ocean contrast and topography 


In order to render our simulations more realistic, we also initialized them with a localized 
wave-packets of the unstable modes. An example of such initial condition is given in Figure 
12. It is obtained by taking a Gaussian envelope in the zonal direction of all fields correspond- 
ing to the most unstable mode, with a width roughly corresponding to two wavelengths of 
the latter. (It is possible, by solving a “dry” adjoint problem, to identify the most rapidly 


growing non-modal perturbation, following e.g. Parker (2008) in the present context, but we 
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Figure 9. Evolution of the most unstable mode in moist-convective environment, as seen in the vorticity and velocity 
fields in the upper (upper row) and lower (second row) layer. Corresponding evolution of condensation (third row) and 
divergence (fourth row) in the lower layer. Same simulation as in Figure 8, but with condensation and evaporation 
switched on (Colour online). 
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Figure 10. First row: Late stages of the evolution of vorticity and velocity in the upper (left panel) and lower (right 
panel) layers illustrating cyclogenesis in the lower layer due to the moist-convection. Second row: Corresponding con- 
densation (left panel) and divergence (right panel) fields in the lower layer. Third row: Corresponding precipitation (left 
panel) and clouds (right panel). To distinguish various details, we made a zoom on two adjacent detaching vortices 
(Colour online). 


adopt here a more straightforward approach). We analyze the evolution of the jet with super- 
imposed perturbation of Figure 12 in “dry” (not shown) and moist-convective environments, 
and determine the role in this process of topography and the land-ocean boundary, which are 
roughly mimicking the West-African plateau, The height of the plateau is taken to be constant 


b = 0.05Hp, and have zonal extension [—0.35 0.35] in non-dimensional terms. The nonlinear 


September 15, 2021 Geophysical and Astrophysical Fluid Dynamics LLEJ’ GAFD'RZ2021°Jun30 


325 


326 


327 


328 


329 


330 


331 


332 


333 


334 


359 


336 


Sof 


338 


20 M. Rostami and V. Zeitlin 


0.08 
0.06 
0.04 


0.02 


Time 
(ep) 





06 -04 -0.2 0 0.2 04 0.6 06 -04 -0.2 0 0.2 04 0.6 
X X 


Figure 11. Hovmdller diagrams of vorticity in the lower layer at y = 0.1 in “dry” (left panel) and moist-convective 
(right panel) environments corresponding to Figs. 8 and 9, respectively. Note that colorbars are not the same in the 
panels, for better visibility (Colour online). 
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Figure 12. Initial condition as a wave packet constructed with the most unstable mode. Left (right) panel: vorticity 
and velocity in the upper (lower) layer. Global configuration is the same as in Figure 4. The maximal non-dimensional 
amplitude of the perturbation is 0.005 (Colour online). 


simulations are initialized with wave-packets of unstable modes over the land, as this is the 
case in the observations of AEW. The evolution of the wave-packet in the moist-convective en- 
vironment in such configuration is presented in Figure 14. ‘The corresponding development of 
clouds and precipitation is shown in Figure 15. As follows from the Figures, the wave-packet, 
as expected, is subject to dispersion, accompanied by wave emission. At the same time, weak 
lower-layer vorticity anomalies induced by the initial condition are growing in intensity, due 


to the instability of the carrier mode, cf. Figure 9. There is practically no difference between 





the evolution with and without topographic elevation in the dry case, except for a weak to- 
pographic Kelvin wave propagating along the eastern coast (the corresponding signal, which 
is enhanced in the presence of moist convection, can be distinguished in the right-bottom 
panel of Figure 14), as the “continental plateau” is low - as already explained, we consider 
a configuration which is similar to the real one, within the precision of the model. However, 
the land-sea contrast in the moist-convective case is important. Cyclonic vorticity anomalies 


due to initial perturbation move over the “continent” practically without condensation, but 
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Figure 13. Hovmdller diagrams of vorticity (left panel) and condensation (right panel) in the lower layer corresponding 

to vorticity anomalies crossing the land-sea boundary (Colour online). 





d 


when they cross over the ocean the condensation due to evaporation switches on and leads 
to their rapid intensification. The result of this process is an appearance of a series of intense 
vortices, detaching from the continent and moving westward. ‘The Hovmoller diagrams of vor- 
ticity and condensation presented in Figure 13 illustrate this process. Due to periodicity in 


zonal direction, these intense vortices eventually hit the continent at its eastern side. In the 





idealized setup of the present experiments, with the easterly jet of uniform intensity extending 


through the whole zonal extension of the domain, the encounter with the eastern continental 





boundary happens before the vortex has time to clearly detach from the jet, as in Figure 10. 





In real life, the intensity of the jet decreases across the Atlantic, and the detached vortices 
propagate North-West on the G—plane. The evolution of the hurricane-like vortices of different 
intensities was recently studied within the mcRSW model (Rostami and Zeitlin 2020), so we 
refer to this paper for the details of the behavior of the detached cyclones at later stages. 
We should emphasize that, in the context of AEJ, the value of the group velocity of the 
perturbation is important, as it defines the speed of propagation of the wave-packet. We 
repeated the simulations of Figure 14 with a wave-packet built with the mode k = 13 with 
the fastest of all unstable modes group velocity, cf. Figure 3. The results of this, and other 
experiments, where we combined different unstable modes to construct initial wavepackets are 


qualitative and quantitatively similar (not shown). 


5. Summary and Discussion 





We investigated development of instabilities of an easterly jet configuration in the presence of 
moist convection and land-sea contrast in a simple atmospheric model, which allows for low- 


cost high-resolution linear and nonlinear numerical analyses. We established linear stability 
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Figure 14. Upper row: Evolution of vorticity and velocity in the presence of “continental plateau” of non-dimensional 
height b = 0.05, and of moist convection in the upper layer (left panel) and lower layer (right panel). Lower row: The 
corresponding evolution of condensation (left panel) and divergence fields (right panel). Simulations were initialized with 
the wave packet of Figure 12, no radiative relaxation. Notice a difference in the colorbar variation of vorticity field. 
Vertical dashed lines represent the borders of the “continent”. Time= 10, 20, 30 from up to down (Colour online). 





diagrams for different parameters of the jet, and identified the unstable modes. We showed 
that the mechanisms of nonlinear saturation of the instability of the jet with respect to the 
perturbation in a form of the most unstable mode are significantly different in adiabatic 
and moist-convective environments. We demonstrated that saturation of the instability in the 
“dry” case is due to subharmonic interactions and, thus is sensitive to boundary conditions. In 
the moist-convective case we demonstrated a robust generation of a secondary westerly jet and 
of localized zones of intense cyclonic vorticity in the lower layer at the northern flank of the jet, 


similarly to observations. It is the structure of the most unstable mode of the jet in the lower 
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Figure 15. Evolution of high levels of precipitable water in the lower layer (left panel), with “continental plateau” of 
non-dimensional height b = 0.05 delimited by dashed lines, and the corresponding evolution of precipitation (right panel) 
(Colour online). 

layer, with its meridional dissymmetry, combined with the effect of enhancement of cyclonic 
vorticity by moist convection which are at the origin of the cyclogenesis process. In a more 
realistic case of initialization of numerical experiments with localized wave-packets of unstable 
modes over topography mimicking a West-African plateau, we demonstrated that intense 


localized vortices are formed in the moist-convective environment over the ocean, to the west 


of the “continent” and move westward. The low flat topography does not significantly influence 





the scenarios of evolution of the perturbations of the jet, except for a weak topographic 
Kelvin-wave signal. On the contrary, the changes of humidity and evaporation at the land-sea 
boundary are important for enhancement of the perturbations. 

Thus, in spite of a highly idealized setup the simple two-layer moist-convective rotating 
shallow water model does capture essential features of the development of African Easterly 
Waves, and allows to get insights on the fundamental dynamical mechanisms underlying the 
observed structures. 
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Figure Al. Stability diagram of the upper-layer easterly Bickley jet on the f-plane with Bu = 5.7 and Bu = 0.7 in 
terms of phase velocity cp (left panel) and growth rate o (right panel) as functions of the zonal wavenumber k. The first 
and the second unstable modes are represented by black and white squares for Bu = 5.7 and by gray and white circles 
for Bu = 0.7, respectively. Ro = 0.0467, Up = 0.03, a = 12, Hoi = Hog = 0.5, 6 = 0.5. 
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Figure A2. Geostrophic stream-functions of the most unstable mode on the f-plane in the upper layer (left panel) 
[we = ni + s7y2] and the lower layer (right panel) [~1 = 71 + 72]. Positive (negative) values are represented by solid 
(dashed) lines, and the velocity field is represented by arrows. k = 20, Ro = 0.15, Uo = 0.045, a = 20, Hoi = 1/3, Ho2 = 





2/3, = 13° N, 6 =0 (Colour online). 


Appendix A: Stability of the easterly jet on the f- plane 


In this Appendix we present the results of the auxiliary investigation of the stability of the 
easterly upper-layer jet on the f- plane. ‘The stability diagram, and the phase-portrait and 
meridional cross-section of the most unstable mode are displayed, respectively in Figures A1, 
A2, and A3. We present in Figure Al the results both for large and small Burger numbers of 
the jet. The former case is close to that of Lambaerts et al. (2012), with the Rossby number 


being close, too, and the results quantitatively and qualitatively similar, which is, thus, a 





successful benchmark. ‘The latter case highlights the dependence of the results on Bu, which 


was not analyzed in Lambaerts et al. (2012). 
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Figure A3. Meridional cross-sections of u;,v;,; of the most unstable mode on the f-plane shown in the Figure A2, in 


the upper layer (left panel) and the lower layer (right panel) (Colour online). 
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